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Abstract 



At the large transverse momentum region, the production of heavy-heavy 



^ ■ bound-states such as charmonium, bottomonium, and be mesons in high en- 



ergy e + e~ and hadronic collisions is dominated by parton fragmentation. We 
calculate the heavy quark fragmentation functions into the D-wave quarko- 
nium and be mesons to leading order in the strong coupling constant and in 
the non-relativistic expansion. In the be meson case, one set of its D-wave 
states is expected to lie below the open flavor threshold. The total fragmenta- 
tion probability for a b antiquark to split into the D-wave be mesons is about 
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2 x 10 5 , which implies that only 2% of the total pseudo-scalar ground state 
B c comes from the cascades of these orbitally excited states. 

I. INTRODUCTION 

The production of heavy quark-antiquark bound states like charmonium, bottomonium, 
and the yet undiscovered charmed beauty mesons at high energy e + e~ and hadronic ma- 
chines can provide very interesting tests of perturbative QCD. In the charmonium case, our 
traditional wisdom has been taken the following scenario: most of the J/ip comes predomi- 
nately from either the radiative decay of the P-wave Xcj states produced by the lowest order 
gluon fusion mechanism or the weak decays of B mesons, whereas the ip' should be produced 
almost entirely from B meson decays. Similarly, most of the T should be produced from the 
radiative Xbj decay. Only until recently can this naive thinking be confronted by the wealth 
of high energy experimental data collected at the Tevatron. During the 1992-1993 run, the 
CDF detector recorded surprisingly large production rates of the prompt J/ip, ip', and Xcj 
which are orders of magnitude above the theoretical lowest order predictions [[[]. In 1993 
0, it was realized that fragmentation can play an important role in quarkonium production 
at the large transverse momentum (p^) region at the Tevatron. Parton fragmentation into 
quarkonium is formally a higher order process but nevertheless it can be enhanced at the 
sufficiently large transverse momentum region compared with the usual gluon fusion process 
0. A simple explanation of the enhancement is that the fragmenting partons are produced 
at high energies but with small invariant masses which can lead to the enhancement factor 
of powers of (pt/™c) 2 compared with the gluon fusion mechanism. 

Most of the fragmentation functions for quarkonium that are relevant to phenomenology 
have now been calculated. Gluon fragmentation functions were calculated in Ref. @ for the 
S-wave, Ref. |||| for the P-wave, and Ref. || for the 1 D2 state; heavy quark fragmenta- 
tion functions were calculated in Ref. |||7j for the S-wave, and Ref. [§dTIJ for the P-wave. 
Phenomenological applications of these fragmentation functions have been performed by a 
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number of groups pl| , [12| . The current scenario of charmonium production at high trans- 
verse momentum is as follows: among all fragmentation contributions that are relevant to 
J/if> production, the dominant ones come from (i) gluon fragmentation into a color-octet 
3 Si cc state followed by a double El transition into J/ip, and (ii) gluon fragmentation into 
Xcj followed by the radiative decay Xcj — ► J/if > + 7- The former mechanism is the most 
dominant one for if)' production because of the absence of Xcj(2P) states below the DD 
threshold. The current experimental status with this new fragmentation insight has now 



been reviewed by many people in various occasions. We refer the readers to Refs. HTTH T4 
for more details. 

The charmed beauty be mesons are expected to be observed soon at the Tevatron and 



preliminary limits on the production rate already exists ||15|| . Parton fragmentation has 
been applied phenomenologically to the production of be mesons as well at the Tevatron 
|T6| , p!7[| . The heavy quark fragmentation functions for b* — > (be) were calculated in Ref. [18 



for the S-wave and in Refs. |8],|9|] for the P-wave. The validity of the parton fragmentation 
has been questioned in some recent exact 0(af) calculations ]B|. However, these calcu- 
lations obtained different results and led to conflicting conclusions among themselves. It 



was only until the recent appearance of the results of Ref. J2[J that this controversy can 
be settled down and the validity of the fragmentation approximation can be established in 
the production of (be) mesons at the Tevatron. One can, therefore, use the fragmentation 
approximation with confidence in calculating the production cross sections and transverse 
momentum spectra of the (be) mesons, including the leading logarithmic corrections, induced 



gluon fragmentation contribution |17j , and contributions from all orbitally excited states. 

The purpose of this paper is to extend the previous calculations of heavy quark frag- 
mentation functions to the D-wave case. The D-wave orbitally excited states are of interests 
phenomenologically. In the (be) case, one set of the D-wave states is predicted by potential 
models to lie below the BD threshold. These states, once produced, will cascade into the 
pseudo-scalar ground state B c via pion and/or photon emissions and thereby contribute 
to the inclusive production rate of B c . The excited D-wave charmonium resonances have 



also been suggested |2l| to resolve the if)' surplus problem observed at CDF. In Section II, 
we briefly review the factorization model of Bodwin, Braaten, and Lepage [^2j that can 
consistently factor out the long-distance physics and short-distance perturbative factors for 
the inclusive production and decay rates of heavy quarkonium. In Section III, we present 
the calculation of the D-wave fragmentation functions for both the unequal and equal mass 
cases within the spirit of the factorization model. We discuss our results in Section IV and 
conclude in Section V. 

II. THE FACTORIZATION MODEL 

A rigorous theory based on non-relativistic QCD (NRQCD) has been developed recently 



by Bodwin, Braaten, and Lepage |22j for the inclusive decay and production of heavy quarko- 
nium. Using this formalism, the fragmentation function for a heavy quark Q to split into a 
quarkonium state X with longitudinal momentum fraction z can be written as 

D Q ^ x (z,fi) = Y,d n (z^)(OX) , (1) 

n 

where O n are local 4-fermion operators defined in NRQCD. The short-distance coefficients 
d n (z,ii) are independent of the quarkonium state X. For a fragmentation scale /i of order 
of the heavy quark mass itlq, the coefficients d n (z,/j,) can be calculated using perturbation 
theory in strong coupling constant a s (2mg). The relative size of the matrix elements (O^) 
for a given state X can be estimated by how they scale with tuq and with the typical relative 
velocity v of the heavy quarks inside the quarkonium. Thus, the factorization formula Eq. ([!]) 
is a double expansion in a s and v. To determine the relative importance of the various terms 
in this formula, one should take into account both the scaling in v of the matrix elements 
(0*) and the order of a s in the coefficients d n (z,fi). The leading term in v of this formula 
corresponds to the popular color-singlet model [p]|1 . However, keeping only the leading term 
in this double expansion can sometimes lead to incomplete or even inconsistent results due 
to the presence of infrared divergences . 
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One important feature in the factorization model is that the quarkonium state X is 
no longer considered as solely a QQ pair but rather a superposition of Fock states. For 
example, the spin-triplet 3 Dj quarkonium states, denoted by 5j, have the following Fock 
state expansion, 

|#) = 0(l)\QQ( 3 Dj,l)) + 0(v)\QQ( 3 Pj,,8)g) + 0(v 2 )\QQ( 3 S 1: 8 or l)gg) + ■■■ (2) 

where the notations 1 and 8 refer, respectively, to the color-singlet and color-octet states 
of the QQ pair. In the above Fock-state expansion there are also other 0(v 2 ) states, e.g., 
IQQ^Dj, 8 or l)gg), but their production will be further suppressed by powers of v. Heavy 
quark fragmentation into D-wave quarkonium Sj can be deduced from the leading Feynman 
diagrams of Q* — > (QQ)Q as shown in Fig.l and Fig.2. In Fig.l, the QQ pair can be either 
in the color-singlet 3 Dj state, in the color-octet 3 Pj', or in the color-octet or color-singlet 
3 Si states, whereas in Fig.2 it can only be in the color-octet 3 Si state. The short- distance 
factors deduced from these diagrams are all of order a 2 . Furthermore, all these diagrams 
are also of the same order in v. According to the NRQCD ||22|| , the relevant local 4-fermion 
operators describing these three Fock states are Oi( 3 Dj), Os( 3 Pj'), and Os( 3 Si) which have 
relative scalings like v 4 , v 2 , and v°, respectively. However, the QQ pair in the color-octet P- 
wave and S-wave states can evolve nonperturbatively into the physical D-wave quarkonium 
state 5 j by emitting one and two soft gluon(s), respectively. Emitting a soft gluon costs a 
factor of v at the amplitude level in NRQCD and hence a factor of v 2 in the probability. 
Thus, with a non-relativistic normalization convention in the state, all three matrix elements 
(Ol J ( 3 Dj)), (0 8 J ( 3 PjO), and (0 8 J ( 3 Si)) scale like m 3 Q v 7 . As a result, both Fig.l and Fig.2 
contribute to the fragmentation functions of order a 2 v 7 and should all be taken in account for 
a consistent calculation. With appropriate normalization, the color-singlet matrix elements 
(Ol J ( 3 Dj)) as well as (Of( 1 D 2 )) for the spin-singlet D-wave state 5® can be related to 



the non-relativistic radial D-wave wave-function Rp(Q), which is the spin average of the 
spin-singlet and spin-triplet i? D (0)'s, by 

(OfCD 2 ))^^\RW)\\ (3) 
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(of?f Cj)> K 15(2^ ^ (4) 

Potential models can be used to determine the value of the wave-function so that the NRQCD 
matrix elements for the color-singlet contributions are fixed. Unfortunately, potential mod- 
els cannot determine the color-octet matrix elements since dynamical gluons are involved. 
However, from our experience of the P-wave quarkonium case ||, we do not expect the 
color-octet component to play a major role in the heavy quark fragmentation. We note that 
this is in sharp contrast with the gluon fragmentation in which a gluon can fragment into 
the color-octet 3 5i state via the process g* — > QQ. This process is of order a s and is at least 
one power of a s lesser than the leading color-singlet term in the gluon fragmentation. We 
also note that the color-singlet contribution of heavy quark fragmentation from Fig.l is free 
of infrared divergences. On the other hand, infrared divergences will show up in the gluon 
fragmentation function into the spin-triplet D-wave quarkonium. It is necessary to include 
the color-octet contributions in order to achieve a finite and sensible perturbative answers 
for the gluon fragmentation functions into spin-triplet D-wave states. In what follows we 
will restrict ourselves to the color-singlet contribution in the heavy quark fragmentation into 
D-wave quarkonium. If the color-octet matrix elements can be determined in the future, 
it is straightforward to include their contributions since their corresponding short- distance 
coefficients can be extracted easily from the previous S-wave and P-wave calculations by 
modifying the color factors. 

The factorization model for the quarkonium system can be extended to the unequal mass 
case like the (be) meson system. In what follows, we will denote the spin-singlet and the 
spin-triplet D-wave (be) meson states by 5 bc and S b f, respectively. In this case, figure 2 is 
absent for the color-octet contribution. 
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III. HEAVY QUARK FRAGMENTATION FUNCTIONS INTO D-WAVE 

HEAVY-HEAVY MESONS 



The general covariant formalism for calculating the production and decay rates of S-wave 
and P-wave heavy quarkonium in the non-relativistic expansion was developed some times 



ago |25]. It is straightforward to extend it to the case of unequal mass and higher orbital 
excitation. We shall present a covariant formalism for the production of D-wave meson 
(Qq) by some unspecified short-distance processes. Let ttiq and m q be the masses of the 
two quarks with tuq > m q , and introduce the mass parameters r = m q /(mQ + m q ) and 
f = mQ/(rriQ + m q ) = 1 — r. In the leading non-relativistic approximation, the mass M 
of the meson is simply tuq + m q . The amplitude for producing the bound-state (Qq) in a 
state with momentum P, total angular momentum J, total orbital angular momentum L, 



and total spin 5* is given by 



r r/ 3 k 

MP) = E J 7^1 ( ll z; SS z \JJ z ) M(P, k) , (5) 



Lz^Sz 

where 



(27T) 



M(P,k) = O r Tss z (P,k). (6) 

Or represents the short- distance interaction producing the Q and q in a specific relative 
orbital angular momentum, and in general is a product of Dirac matrices. The heavy quarks 
Q and q have momentum rP + k and fP — k, respectively, where 2k is the small relative 
momentum of the quarks inside the bound state, ^^(k) is the Bethe-Salpeter wave- 
function in the momentum space and is assumed to be a slow-varying smooth function of k. 
Sometimes A4(P, k) is a trace by which the on-shell spinors v(rP + k, s) and u(fP — k, s) 
are brought together in the right order. Tss z (Pi k), up to second order in k, is given by 



r SSz (p,k) = J mQ + m « WVi 



/ ( rP + k — m 

V W 3 + r M 2m q 



SSz } v(rP + k,s) u (rP — k, s) , 
ff- k' + rriQ^ 



( ,5 \ 



\-t{P,S z ) J 



2m Q 



(7) 



where in the middle parenthesis 7 5 is for 5 = Sz = (spin-singlet) and —^(P,Sz) is for 
5 = 1 (spin-triplet). 

Since k/M is a small quantity, we can expand A4(P,k) around k = in a Taylor 
expansion: 

dM(P,k) 



M{P,k) =M(P,0) + k o 



dk n 



k=0 



1 dlM(P,fc) 
+ 2 Q ^ dk a k p 



+ 



(8) 



k=0 



where the first, second, and third terms correspond to quantum numbers L — 0, 1, and 2 of 
the orbital angular momentum, and so forth. Thus for the S-wave (L = 0), P-wave (L = 1), 
and D-wave (L = 2) states, the amplitude A(P) will depend on the radial wave-functions 
-Rs(O), -Rp(O), and -R^(O) through the following relations: 

d3k *„(k) = ^ (0) 



(2tt 

— -p ^(k)^ = -i^—R P {Q)e a {P,L z ) 

<i 3 k / 15 

^(k)^^ = \ — R' D (Q)e a p(P, L z ) , 



(2tt 



(9) 



where e Q is the polarization vector for the spin-1 particle, and e Q/ g is the totally symmetric, 
traceless, and transverse second rank polarization tensor for the spin-2 particle. In the 
D-wave case, the amplitude becomes 



2 V 07T OkaOkp 



(10) 



fc=0 



for the spin-singlet case, where Jz = L z , and 

1 /l5 



(11) 



fc=0 



where 



n^ p (P,J z )= J2 e af3 (P,L z ) e p (P,Sz)(2L z ;lS z \JJz) 

Lz,Sz 



(12) 



for the spin-triplet case. Using the appropriate Clebsch-Gordan coefficients, we have 
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Kf P {P, J z) = ~\J^ (I'Pap e„(P, Jz) - V ap € P (P, Jz) - Pp P e a (P, J z )) , (13) 

Kf P (P Jz) = ^= (e aa (P, J z )e TPpa , P T g™' + e Pff (P, J z )e Tap(T , P T g™') , (14) 

n J a f p (P,J z ) = e a p p (P,J z ), (15) 

where 

Va(3 = -dap+^f , (16) 

and e a/ 3 P (P, J z ) is the totally symmetric, traceless, and transverse spin-3 polarization tensor. 
The polarization sums for J = 1, 2, and 3 are given by the following familiar expressions 



£ e a (P,J z )e* fj (P,J z )=V a p, (17) 

Jz=-1 

2 1 1 

E ^( P > ^)<>(^ Jz) = g (^V/to + ^a^p) " g^T'p.r , (18) 



Jz=-2 
3 



Jz=-3 



2J e af3-y(P, Jz)£* parj {.Pi Jz) ~ 7T ( PapPpaP^-q + PapP^P-ya + P otoP j3 pP -frj 



~^~P aaP fir\P ^ p ParjPfjaP-yp Par/PfipPja 

i^P affP ^r)P per ~\~ PaftP^ycrPpr) PafsP^ypPar/ 

~\~P orfP PrfP pa PayPf3<rPpri 'aj'fip'arj 

^P ft~{P ocqP pa P (}~{P aaP pr] P(3"/PapPar^j • (1^) 

After writing down the formalism for the general production of D-wave mesons, we shall 
specify the production mechanism. The specific production of D-wave (Qq) mesons that we 
are considering is by the fragmentation of a heavy quark Q: 

Q*(q)^Qq(P,l) + q(p'), (20) 

where the off-shell Q* is produced by a high energy source V. The fragmentation function 

is then given by the formula |||l^] : 

1 r , „ / M 2 r 2 M 2 \ , E\A(P)\ 2 , s 

D(z) = ds6[s lim ^ . \ ' 1 , (21) 



16n 2 J \ z 1 — z I ^ L4 Q 



where S l^ol 2 = ^cTr(rr^) is the tree-level amplitude squared to create an on-shell Q 
quark with the same 3-momentum q, and J2 \A(P)\ 2 is the amplitude squared for producing 
(Qq) + q from the same source T. The general procedures to extract the fragmentation 
function have been described in detail in Ref. PJTBf and we shall be brief in what follows. 



A. Unequal Mass Case 

We shall first derive the fragmentation functions for a heavy quark Q into an unequal 
mass D-wave (Qq) meson. The fragmentation of Q* into a Qq meson is given by the process 
in Eq. (pO"l), of which the leading Feynman diagram is given in Fig. |||. In reality, the above 
fragmentation process applies to b* — > 5 bc + c (b* — ► 5 b f + c) and also its charge conjugate 
b* — > 5 bc + c (b* — > 5j c + c). The amplitudes that are presented in the following are for 
Q* ~ > (Qq) + q, while the amplitudes for the conjugate process can be obtained by complex 
conjugation. However, the final expressions for the fragmentation functions are valid for 
both Q* — > (Qq) + q and Q* — > (Qq) + q fragmentation processes. 

The amplitude A(P) for producing the spin-singlet D-wave state is given in Eq. fllQ"|) 



with 



d 2 M 



dk a dkp 



given by 

k=0 



d 2 M 



dkadkp 



= N l3 A F u(p')V a ^ 5 T , (22) 

k=0 



where 



N ^ Ac ^^km w,th Cr = l' (23) 

yap = v o0 + (i_ 2r )vf , (24) 



and 



yf = ^-MA%q a q P (P + 2M)(i - fM) + ^-MA 2 F A F q a q? d (P- M) 

2 Af 
+ —^ApApifi (gV + gV)f + —MA F A 2 F (q a n f3 + q p n a ) i (P - M) 

+ -A 2 F i/i (n a ^ + n^ a )P+ 4fMA F n a n 13 t/i iff— M) , (25) 
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Vf = -^MA|( 7 V + W)(4 - fM) - ^MA F A F ^ a qP + 7 V) 

- -MAH^n 13 + ^n a ) , (26) 



with 



Ap = (s — f'M 2 )" 1 , (27) 
A F = (n • (g - fP)) _1 . (28) 

In the above formulas, s = q 2 and we have chosen the axial gauge associated with the 4- 
vector n M = (1, 0, 0, —1). In this gauge, figure is the only leading Feynman diagram where 
factorization between the short distance process and fragmentation becomes manifest! We 
square the amplitude A(P) and sum over colors and helicities of q and (Qq), and put them 
into Eq. ([H]). The heavy quark fragmentation function for the l D 2 state thus obtained is 
given by 

/N a 2 (2rM)\R" D (0)\ 2 1 z(l - zf { , 2 , 

+ 5(267 - 600r + 580r 2 - 160r 3 + 64r 4 )^ 2 

+ 10(-255 + 825r - 1130r 2 + 596r 3 - 88r 4 - 64r 5 )^ 3 

+ (3225 - 13050r + 22585r 2 - 16872r 3 + 5036r 4 - 576r 5 + 1152r 6 )/ 

- 4f(690 - 2535r + 4200r 2 - 2279r 3 + 968r 4 - 524r 5 - 160r 6 )z 5 
+ f 2 (1545 - 4890r + 8525r 2 - 2176r 3 + 3076r 4 - 320r 5 + 320r 6 )/ 

- 2f 3 (255 - 660r + 1475r 2 + 174r 3 + 740r 4 + 120r 5 )z 7 

+ 5f 4 (15 - 30r + 107r 2 + 56r 3 + 80r 4 )/| . (29) 

For the spin-triplet case, the amplitude A(P) is given in Eq. ([TTD with °' M ' 



d 2 M p {k) 



given by 

k=0 



dk a dkp 

where 



= NijApuip^V^r , (30) 

k=0 



yaPp = yapP + ^ _ 2r )Vf P , (31) 
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and 



af3p 



-A F 

r 



(Yg pa + l a 9 P n - -A F ( (q a (YY - 2fg pP ) + q p (YY - 2fg f 



+2M (q a g pl3 + q p g pa ) ) + —M 2 A 2 F q a q (3 Y 



+ -A F i 2{Yg pp + tV^ 
r 



a + fM) 



+A F (n a {{YY - 2rg pl3 )f+ 2rMg pp ) + n p {{YY ~ 2rg pa )f + 2rMg pa ) 



+-A F [q a ((YY - 2rg PIJ )f+ 2rMg pf) ) + q"((YY ~ 2rg pa )f+ 2rMg 



pa • 



AfMA F 



- q a q p A 2 F + —A F A F {q a n p + q p n a ) + A 



2 n a n? 



V, 



a(3p 



M 



Ap ifi 



Ap (q a Y + <?V) + A F (n a Y + n V) 



Y 



M)Y , (32) 
(33) 



We then square the amplitude and sum over the colors and helicities of q and (Qq) using the 
polarization sum formulas in Eqs. flTTD - (0). The fragmentation functions thus obtained 
are 

a 2 (2rM)|i^(0)| 2 i z{ i- z f 



Dq^q^dY^ 3240ttM 7 r 6 f 2 (1 - fz) w 

x |60(1 + 4r) 2 - 60(1 + 4r)(7 + 36r - 28r 2 )^ 
+ 5(261 + 2652r + 6128r 2 - 12576r 3 + 7360r 4 )^ 2 
+ 10(-237 - 2589r - 8704r 2 + 27328r 3 - 28608r 4 + 11360r 5 )z 3 
+ (2775 + 31350r + 153935r 2 - 637872r 3 + 897816r 4 - 611056r 5 + 183552r 6 )z 4 

- 4f(540 + 6585r + 48660r 2 - 162629r 3 + 199928r 4 - 118084r 5 + 28400r 6 )z 5 

+ f 2 (1095 + 13830r + 137935r 2 - 376916r 3 + 385696r 4 - 192320 5 + 36800r 6 )z 6 

- 2f 3 (165 + 2100r + 26395r 2 - 55836r 3 + 48360r 4 - 26720r 5 + 3360r 6 )z 7 

+ 5f 4 (9 + 114r + 1717r 2 - 2424r 3 + 2476r 4 - 1872r 5 + 192r 6 )/ J , (34) 



, , a 2 (2rM)\R" D (0)\ 2 1 z(l - z) 2 . 
Dq^d 2) (z) = A a J:J 7 K n ^ J m \ 180 + 180(-7 + Ar)z 



648vrM 7 r 6 f 2 (1 - fz) 10 
+ 5(777 - 852r + 464r 2 - 64r 3 + 128r 4 )z 2 

+ 10(-693 + 1087r - 1126r 2 + 624r 3 - 112r 4 - 128r 5 )^ 3 
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+ (7875 - 15650 + 23295r 2 - 22564r 3 + 9112r 4 + 128r 5 + 2304r 6 )/ 

- 4f(1470 - 2005r + 4740r 2 - 4613r 3 + 2296r 4 - 488r 5 - 320r 6 )^ 5 
+ f 2 (2835 - 2050r + 12035r 2 - 8332r 3 + 5952r 4 - 960r 5 + 640r 6 )z 6 

- 2f 3 (405 - 40r + 2565r 2 - 602r 3 + lOOOr 4 - 320r 5 )^ 7 

+ 5f 4 (21 + lOr + 205r 2 + 52r 3 + 84r 4 )/| , (35) 
2a 2 (2rM)|^(0)| 2 1 z{\ - zf f 

+ 5(399 - 174r + 200r 2 )^ 2 + 10(-378 + 193r - 437r 2 + 70r 3 )z 3 

+ (4725 - 2900r + 8390r 2 - 1918r 3 + 2268r 4 )^ 4 

+ 2(-1995 + 1700r - 4930r 2 + 1617r 3 - 2387r 4 + 350r 5 )^ 5 

+ (2205 - 2830r + 7940r 2 - 5194r 3 + 4914r 4 - 930r 5 + lOOOr 6 )/ 

- 2f(360 - 325r + 1680r 2 - 749r 3 + 1540r 4 + 295r 5 + 90r 6 )/ 

+ 5f 2 (21 - 14r + 133r 2 - 56r 3 + 189r 4 - 18r 5 + 18r 6 )z 8 } . (36) 

All the above fragmentation functions are valid for the fragmentation of b — ■> 5 bc and 5 b f, as 
well as for b -> W c and 5f, where J = 1, 2, 3. 

For a given principal quantum number n the 1 D 2 and 3 D 2 states constructed in the LS 
coupling scheme are mixed in general to form the physical states 2 + and 2 + ' defined as 

|2+') = cos9\ 1 D 2 ) -sin#| 3 £> 2 ) , 

(37) 

|2+) = sin9\ 1 D 2 ) + cos9\ 3 D 2 ) , 

where 6 is the mixing angle. Thus, the fragmentation functions into the physical states 2 + 
and 2 + ' are 

D Q^Qq(2+')( z ) = cos 2 eD Q ^ Q g ( i D2) (z) +sm 2 6D Q ^ Q q { t D2) {z) - sin cos 6 D mix {z) , 

(38) 

D Q ^ Q q( 2 +)(z) = sin 2 QD Q ^ Qqi i D2) (z) + cos 2 9D Q ^ Qq - ( 3 D2) (z) + sin 9 cos 9 D mix (z) . 

Therefore, we also need to calculate the fragmentation function D mix (z). We note that in the 
case of charmonium it is the charge-conjugation quantum number C that prevents the mixing 
of cc( l D 2 ) and cc( 3 D 2 ) states, because cc( 1 D 2 ) has C — +1 while cc( 3 D 2 ) has C = —1. In 

13 



the unequal-mass case the meson does not have this quantum number C to prevent the 
mixing. The mixing fragmentation function is obtained by calculating the interference term 
A( 1 D2)A*( 3 D 2 ) + A*( 1 D2)A( 3 D2) and sum over the colors and helicities of the particles in 
the final state as before. The mixing fragmentation function thus obtained is 



Anix(^) = 



5a 2 s (2rM)\K 



LA "" 2 1 : " lf 1 12 ■ L2(-5 • 4r- 12,-): 



5AV6nM 7 r 6 f 2 (l-fzf 
+ (129 - 198r + 156r 2 + 64r 3 )^ 2 + 4(-39 + 85r - 76r 2 + 8r 3 - 36r A )z 3 

+ 2(57 - 156r + 133r 2 - 130r 3 + 198r 4 + 48r 5 )/ 

- 4f(12 - 27r + 4r 2 - 51r 3 - 2r 4 - Ur 5 )z 5 

+ f 2 (9 - 16r - llr 2 - 62r 3 - 4r 4 )/| . 



(39) 



The fragmentation probabilities for these states can be obtained by integrating the frag- 
mentation functions over z, 

ION 



D 



Q^Qq{ 1 D 2 ) 



[z)dz 



f(793 + 3712r + 42160r 2 + 123220 3 + 253315r 4 



2i r 7 f io 

+ 126622r 5 + 9778r 6 + 1180r 7 - 80r 8 ) 

+ 315r(ll + 40r + 280r 2 + 580r 3 + 701r 4 + 152r 5 + 16r 6 ) logr 



(40) 



r 1 N 
J o D Q ^Q 9{3Dl) (z)dz = 



J D Q-Q^D 2) {z)dz=^-^ 



f(209 + 457r + 28224r 2 + 479388r 3 + 1136191r 4 

- 358043r 5 - 430066r 6 - 113308r 7 + 12528r 8 ) 

+ 105r(5 + 60r + 1772r 2 + 8644r 3 + 4759r 4 - 5632r 5 - 2108r 6 

- 368r 7 + 64r 8 )logr 
5N 



(41) 



1127 - 700r + 67880r 2 + 320552r 3 + 393317r 4 



+ 191990r 5 + 118742r 6 + 4712r 7 - 160r 8 ) 

+ 315r(9 + 44r + 604r 2 + 1260r 3 + 879r 4 + 532r 5 + 156r 6 ) logi 



(42) 



(z)dz 



IN 



'2191 + 4828r + 100786r 2 + 160027r 3 + 288799r 4 
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+ 428008r 5 - 63854r 6 + 87913r 7 + 14422r 8 ) 

+ 630r(15 + 20r + 343r 2 + 126r 3 + 861r 4 + 112r 5 + 63r 6 

+ 78r 7 + 6r 8 ) logr , 



(43) 



o 



D mix (z)dz 



f(153 + 358r + 3214r 2 - 13796r 3 - 15021r 4 



+ 7176r 5 + 9696r 6 + 240r 7 ) 

+ 105r(5 + lOr - 20r 2 - 180r 3 - 37r 4 + 110r 5 + 36r 6 ) logr 



(44) 



where N = a 2 (2rM)\R' D (0)\ 2 /(3240ttM 7 ). We note that the running scale used in the strong 
coupling constant a s is set at 2rM, which is the minimal virtuality of the exchanged gluon 



B. Equal Mass Case 

Heavy quark fragmentation functions for the D-wave quarkonium can be obtained simply 
by setting r = 1/2 and 9 = in the previous formulas for the unequal mass case. For 
convenience we also present their explicit formulas in the following. 

,,2/n — M d" m\\2 „n „\2 



= 8al(2rn Q )\R (O^zil-zy , _ ^ + _ 
8\ixm n (2 — z) w v 



+ 35328/ - 29344/ + 18344/ - 5848/ + 775/J , (45) 

, N 8a 2 (2m )\R" n (0) I 2 z(l - z) 2 , 9 , 

D n ,q(z) = — ^ Qn 7 DK n / r^r (17280 - 103680z + 321120/ - 551840/ 

Q-V v ; 4057rm^ (2 - 2) 10 V 

+ 546744/ - 314752/ + 112402/ - 24594/ + 2915/) , (46) 

= 16^ 2 (2m Q )|^(0)| 2 z(l - ,) 2 , _ _ 
; 81vrm^ (2-/) 10 V 

+ 58604/ - 38372/ + 16517/ - 4014/ + 445/) , (47) 



/N 8a 2 2m )LR n 2 2 1 -z 2 , 9 , 

D n xQ (z) = — ^ g;i 7 DV ;| / r^r (11520 - 69120z + 231680/ - 488960/ 

+ 675136/ - 592288/ + 309688/ - 80736/ + 8285/) . (48) 



We can also obtain the fragmentation probabilities by integrating over z: 
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f 1 r , s 8a 2 J2m o )\R"n(0)\ 2 ( 776677 \ , s 

/ dzD^g z = sl 7 DK n -^r— - 53355 log 2 , 49 

Jo 8l7rmQ V 21 J 

/ ^^M? (*) = 40 W ( 547127 - 789300 lo S 2 ) > (50) 

Q 

, r, / \ 16^(2/710)1^(0) | 2 Z2889265 n . fl1n1 \ 

I ^<?M = 8 A (^6^ - 24810 '° g2 ) ' (51) 

/' izDq jo(2 , = ^(yi^WI' f ™ _ 3560251o g2 ) . (52) 
Jo Q ^ 5 3 405vrm^ V 21 6 / 



IV. DISCUSSIONS 



We have seen that the covariant formalism used to calculate the D-wave fragmentation 
functions is quite cumbersome and tedious. The complicated results are therefore needed to 
be cross-checked by other means. To be sure we have recalculated the D-wave fragmentation 
functions for the equal mass case from scratch and we did end up with the same results as 
what we have in the above for r = 1/2. Two additional nontrivial checks are performed by 
using the Braaten-Levin spin-counting rule and the heavy-quark spin-counting rule which 
we will discuss in turns. 



A. Braaten-Levin Spin-counting Rule 

The fragmentation function Di_,H(z) for a parton i splitting into a hadron H is related to 
the distribution function fi/n{x) of finding the parton i inside the hadron H by the analytic 
continuation [E71E8 



fi/n(x) = xD^H [-) ■ (53) 



x 



The results of perturbative fragmentation functions allow us to study the perturbative tail 
of the distribution functions of the heavy quark inside the heavy mesons as well. From our 
explicit calculations, we see that fi/H{x) has a pole located at x — f. The pole is cut off 
by non-perturbative effects related to the formation of bound states of the (Qq) pair. This 
pole is of order 6, 8 and 10 for S-, P-, and D-wave states, respectively. In the general case of 
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L-waves we expect this pole is of order 6 + 2L. Therefore, we can expand f(x) as a Laurent 
series, 

tl \ fl "( r ) flr*-l(0 , ■ , fr A \ 

fix) = — + : — - + less singular terms , (54) 

v ' (x-f) n (x- f)"- 1 v ' 

with n = 6 + 2L for the general L-waves. 

The Braaten-Levin rule states that the leading r-dependent coefficients a n (r) satisfy the 

simple spin-counting. Specifically, we have 

a 6( 1 >S'o) : Q6( 3 Si) = 1:3 for S-waves, (55) 

a 8 ( x Pi) : a 8 ( 3 P ) : a 8 ( 3 Pi) : a 8 ( 3 P 2 ) =3:1:3:5 for P-waves, (56) 

a w { x D 2 ) : a 10 ( 3 Pi) : a 10 { 3 D 2 ) : a 10 ( 3 P 3 ) = 5 : 3 : 5 : 7 for D-waves, (57) 
and so on. 

The applicability of these counting rules has been demonstrated for the S- and P-wave 
cases in Ref. [|J. For the D-wave case we can expand the distribution functions, which are 
obtained from the D-wave fragmentation functions by Eq. ([53]), around x = f and obtain 
the first two terms in the Laurent series: 

. 3072(W(rf) 2 320iV(576r - 275)rf 
fQ/M*) = {x _ f yo + r— - 59 + • • • ' ( 58 ) 

. . 18432A^(rf) 2 192iV(576r - 275)rf . . 

f Q /s?^) = {x _ f yo + 7— ~f + " • • • (59) 

, x 30720A^(rf) 2 320iV(576r - 275)rf 

. s 43008A^(rf) 2 448A^(576r - 275)rf , . 



3 



where = a 2 |P^(0)| 2 /(32407rM 7 ). Therefore, the leading coefficients aio's indeed satisfy 
the spin-counting ratio 5:3:5:7. Actually, we found that not only the first term but also 
the second term in the Laurent series of f(x) obey the Braaten-Levin spin-counting rules, 
while the third term does not. We conjecture that the Braaten-Levin spin-counting rule can 
be applied to the first two coefficients in the expansion fl51| ) for all L-wave (L = 0, 1, 2, ...) 
cases. 
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B. Heavy-Light Limit 



Hadrons containing a single heavy quark exhibit heavy quark symmetry in the limit 
rriQ /Aqcd — * oo. In the limit of uiq/Aqcd — > oo, both the heavy quark spin Sq and the 
total spin J of a heavy hadron containing a single heavy quark Q become good quantum 
numbers. This implies that in the spectroscopy of the hadron containing a single heavy 
quark Q, the angular momentum of the light degrees of freedom Ji = J — Sq is also a 
good quantum number. We refer collectively to all the degrees of freedom in the heavy-light 
hadron other than the heavy quark as the light degrees of freedom. For heavy-light (Qq) 
mesons, J\ = S q + L where S q is the spin of the light quark q and L is the orbital angular 
momentum. Thus, hadronic states can be labeled simultaneously by the eigenvalues j and ji 
of the total spin J and the angular momentum of the light degrees of freedom Ji, respectively. 



In general f29fl , the spectrum of hadrons containing a single heavy quark has, for each ji, a 
degenerate doublet with total spins j + = ji + 1/2 and j_ = ji — 1/2. (For the case of ji = 0, 
the total spin must be 1/2.) For D-wave heavy-light mesons, ji can either be 3/2 or 5/2. 
Thus j + ) = (1, 2) and (2,3) for ji = 3/2 and 5/2, respectively. As a result, we expect to 
have two distinct doublets ( 3 -Di, 2 + ') and (2 + , 3 D 3 ) in the limit of m,Qjm q — >• oo, i.e. r — » 0. 
In this limit, the mixing coefficients in Eq. (|3~?D can be determined by the Clebsch-Gordan 
coefficients in the tensor product of a spin 1/2 state and a spin 2 state with the following 
result: 

|2 + ') = yfl 1 ^) + \/f | 3 ^> , (62) 

|2+> = -/f |^ 2 > + y| 3 ^ 2 > , (63) 

i.e., we are transforming the states l D 2 and 3 D 2 in the LS coupling scheme to the states 2 + 
and 2 + ' in the jj coupling scheme. 

In their discussions of the heavy quark fragmentation functions within the context of 



Heavy Quark Effective Theory, Jaffe and Randall PDA showed that fragmentation functions 



can have a I/tuq expansion if expanded in terms of a more natural variable 



rather than the usual fragmentation variable z, and the heavy-quark mass expansion is given 
as a power series in r, 

D(y) = ^a(y) + b(y) + 0(r) , (65) 

where a(y),b(y), etc. are functions of the variable y. The leading term a(y) is constrained 
by the heavy quark spin-flavor symmetry while all the higher order terms contain spin- 
flavor symmetry breaking effects. One can recast our results for the D-wave fragmentation 
functions derived in Sec. Ill in the above form, by carefully expanding the powers of r and 
(1 — fz). The leading 1/r terms of the fragmentation functions are then given by 

10N'(v — I) 2 

D Q ^ S Q,(y) -> — fro ( 3072 - 265Q y + 205Q y 2 - 912 y 3 + 66 °y 4 - l2 °y 5 + 75 ^ 6 ) > 

3N'(v - I) 2 

D Q ^s?<(y) — 10 ( 6144 - 5312 y + 1952 y 2 - 182 V + 102 °y 4 - 12 °y 5 + 15 ^ 6 ) > 

5N'(v — I) 2 

D Q ^ s Q,(y) -> ; (6144 - 5312y + 3392y 2 - 1824y 3 + 1220y 4 - 200y 5 + 105y 6 ), 

1 1 7N'(v — I) 2 

Dq^Mv) ~* To < 384 - m2, J + M V - 114 » 3 + 95 'f 4 - 2 °f 5 + 15 » 6 ) • 



Z,N'(v — \) 2 

d q^2+'(v) - — ^ib-^- ( 6144 - 531 % + 195 % 2 - 182 V + 1020 ^ 4 - 12 % 5 + !% 6 ) > 

gQN'(v — I) 2 

where iV' = a2|i^(0)|7(32407r(rM) 7 ). Thus, at leading-order of l/m Q , we obtain the 
following spin-counting ratios 

D Q^8?*(y) 3 



D Q ^ 2+ '(y) 5 ' 
D Q ^ 2 +(y) _^ 5 

D Q ^(y) 7 ' 

as expected from heavy quark spin symmetry. 



(66) 
(67) 
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V. CONCLUSIONS 



The detections of the D-wave orbitally excited states of heavy-heavy mesons are much 
more difficult than their ground states. So far, none of D-wave quarkonium has been iden- 
tified. The basic reasons are (i) the production rates are small, (ii) small branching ratios 
in subsequent decays, and (iii) the very small efficiencies in subsequent levels of identifi- 
cation. From the potential model calculations, D-wave charmonium are likely to be above 
the DD threshold and hence they will decay predominantly into a pair of D mesons. In 
principle, D-wave quarkonium can be identified partially by their pure leptonic decay (e.g. 
3 Di — ► or by their subsequent decays into the lower-lying states. It is because of the 

subsequent levels of decays that reduce substantially the branching ratios and efficiencies in 
identifying these D-wave bound states. However, both the branching ratios and efficiencies 
are uncertain at this stage. In the following, we shall estimate the production rates of the D- 
wave charmonium, bottomonium, and (be) mesons based upon the fragmentation functions 
obtained in Sec. Ill and IV. 

The values of \R'^(0)\ 2 for the (cc), (bb), and (be) systems were summarized nicely in Ref. 
pl| using different potential models. We choose the QCD-motivated potential of Buchmuller 
and Tye |32[, and the corresponding values of |i?^(0)| 2 for the (cc), (bb), and (be) systems are 
0.015, 0.637, and 0.055 GeV 7 , respectively. These values are for the first set of the D-wave 
states. We note that the values of |i?^(0)| 2 obtained using the Buchmuller- Tye potential 
are the smallest ones of the four potential models studied in Ref. [[31 1. 

In the following analysis, we will adopt the values m c = 1.5 GeV and = 4.9 GeV 



that were used in the fit of the Buchmuller- Tye potential to the bound state spectra 31 



We first evaluate the fragmentation probabilities for a b antiquark into D-wave (be) mesons 
by substituting r = m c /(m b + m c ) = 0.234 into Eqs. (|40|) - ([43]). Using a s (2m c ) « 0.253, 
|i^(0)| 2 = 0.055GeV 7 for the (be) system, and ignoring the possibly small mixing effects, 
we obtain 

/' 

Jo 
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D- b ^ &bc (z)dz = 6.7 x 10~ 6 , (68) 



/ D b ^ &b c{z)dz = 1.7 x 10~ 6 , (69) 

J 

f 1 D- b ^ &bc (z)dz = 6.5 x 1(T 6 , (70) 
Jo 2 

/ D- b ^ S bc(z)dz = 8.5 x 10~ 6 . (71) 

J 

Therefore, the total probability is about 2.3 x 10~ 5 . We can compare this value with the 
corresponding probabilities of the S-wave and P-wave states. Using the same a s and the 
values of |.Rs(0)| 2 and \R' P (0)\ 2 calculated with the same Buchmiiller-Tye potential, the total 
fragmentation probabilities of the IS and IP (be) states are about 1 x 10~ 3 and 1.7 x 10~ 4 , 
respectively |J . These total fragmentation probabilities show a similar suppression factor of 
order v 2 in going from IS to IP states and from IP to ID states, which are in accord with 
the velocity counting rules in NRQCD 

Next, for the D-wave charmonium we use a s (2m c ) » 0.253 and |i^(0)| 2 = 0.015GeV 7 , 
and obtain 

f 1 D c ^ s Jz)dz = 3.1 x 10" 6 , (72) 
Jo 

[ D c ^ 5 c(z)dz = 2.3 x 10~ 6 , (73) 
Jo 1 

D c ^ S c(z)dz = 3.6 x 10~ 6 , (74) 
D c ^ S c(z)dz = 1.7 x 10~ 6 . (75) 

The total probability is about 1.1 x 10~ 5 , to be compared with the probabilities of P(c — > 
J /if)) ps 1.8 x 10 -4 || and P(c — > h c , Xc) ~ 8 x 10~ 5 || obtained using the same inputs. With 
our inputs, the gluon fragmentation probability into the 1 D 2 charmonium || is about 1.2 x 
10 -6 . Similarly, using a s (2m b ) ps 0.174 and \R'^(0)\ 2 = 0.637GeV 7 for the bottomonium, we 



obtain 



D b ^ sb (z)dz = 1.6 x 10~ 8 , (76) 

f 1 D b ^ sb (z)dz = 1.2 x 10" 8 , (77) 
Jo 1 

/ D b ^ 5b (z)dz = 1.8 x 10" 8 , (78) 
jo 2 

D 6 ^(z)(i2 = 8.5xl0- 9 . (79) 
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The total probability for the D-wave bottomonium is very small ~ 5.4 x 10~ 8 . This is to be 
compared with the probabilities of P(b —>■ T) « 2 x [[J and P(b h b , Xb) ~ 2.5 x 10~ 6 . 
We observe that the relative probabilities for the S-wave to P-wave quarkonium and for the 
P-wave to D-wave quarkonium are roughly suppressed by a similar factor of v 2 as in the be 
case. 

At the Tevatron the 6-quark production cross section is of order 10 fib with ?t 6 
GeV, which implies about 10 9 6-quarks with an accumulated luminosity of 100 pb~ x . Using 
the fragmentation probabilities calculated above, we expect about 10 4 D-wave (be) mesons, 
while there are only about 50 D-wave bottomonium. At large Pt, charm quark production 
is very similar to bottom quark production. Therefore, the c-quark production cross section 
is also of order 10 fib and there should be about 10 4 D-wave charmonium for the same 
luminosity. Although the production rates of D-wave (be) mesons and charmonium are 
rather substantial, the small branching ratios and efficiencies of identification in subsequent 
levels of decays render their detections extremely difficult. For example, the combined 
branching ratio Br( 1 D 2 (cc) — > 1 Pi r y, X P\ — > J/ipTc°, J /if} — > fi + fi~) has been estimated to be 
about 10~ 4 []5J. Even after the installation of the Main Injector when the yearly luminosity 
can be boosted to 1 - 3 fb^ 1 , the chance of detecting the D-wave (be) mesons or charmonium 
is still quite slim. The future Large Hadron Collider (LHC) can produce 10 12 — 10 13 bb pairs 
and cc pairs per year running. It can then produce about 10 7 — 10 8 D-wave (be) mesons and 
charmonium. After taking into account the branching ratios and efficiencies, there should 
be a handful number of these D-wave states that can be identified. Nevertheless, detection 
of the D-wave bottomonium does not seem to be feasible even at the LHC! The production 
of bb and cc pairs at LEP, LEPII, and other higher energy e + e~ or fi + fi~ colliders are much 
smaller than at the hadronic colliders. Therefore, it is much more difficult to identify the 
D-wave quarkonium or (be) mesons at the e + e~ or fi + fi~ facilities. 

In conclusions, we have computed the heavy quark fragmentation functions into D-wave 
heavy-heavy mesons containing two heavy quarks to leading order in strong coupling con- 
stant and to leading order in the non-relativistic expansion within the framework of fac- 
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torization model of Bodwin, Braaten, and Lepage. The color-singlet contributions can be 
expressed in terms of the NRQCD matrix elements which are related to the second deriva- 
tive of the non-relativistic radial D-wave wave-function. The color-octet contributions are 
expected to be quite small in the heavy quark fragmentation and therefore we have ignored 
them in the present analysis. The fragmentation functions obtained in this work should be 
useful in the production of the D-wave charmonium, bottomonium, and (be) mesons at the 
future Large Hadronic Collider. 
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FIGURES CAPTION 

1. Feynman diagram for Q*(q) — > QQ{P, 1 or 8) + Q(p'). 

2. Another Feynman diagram for Q*(q) — > QQ(P, 8) + Q(p'). 

3. Feynman diagram for Q*(q) — > Qq(P,V) + q(p'). 
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Fig. 3 



